Macroscopic Quantum Tunneling and Quasiparticle Dissipation in d-wave 
Superconductor Josephson Junctions 
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We examine the macroscopic quantum tunneling (MQT) in high-Tc superconductor Josephson 
junctions with a d-wave order parameter. Using microscopic Hamiltonian and the functional integral 
method, we analytically obtain the MQT rate (the inverse lifetime of the metastable state) for the 
c-axis twist Josephson junctions. In the case of the zero twist angle, the system shows the super- 
Ohmic dissipation due to the presence of the nodal quasiparticle tunneling. Therefore, the MQT 
rate is strongly suppressed in compared with the finite twist angle cases. 
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In the current biased Josephson junctions, the states 
of non-zero supercurrent are metastable owing to tran- 
sitions to lower-lying minima of the potential. At suffi- 
ciently low temperatures, such transitions can be caused 
by macroscopic quantum tunneling (MQT)i*2 through 
the potential barrier. The possibility of observing the 
MQT in Josephson junctions was first pointed out by 
Ivanchenko and Zi'lberman^ The first clear experimen- 
tal observations of the MQT were made in 1981 on small 
s-wave Josephson junctions by Voss and Webb (Nb)£ and 
Jackel et al (Pb)£ 

Since macroscopic systems are inherently dissipative, 
there arises a fundamental question of what is the effect 
of dissipation on the MQT. This issue was first solved by 
Caldeira and Leggett in 1981 by using the path-integral 
method and they showed that the MQT is depressed by 
dissipation. 6,7 This effect has been verified in experiments 
on s-wave Josephson junctions shunted by an Ohmic nor- 
mal resistance Rs (< Rn- the tunnel resistance of the 
junction) As was mentioned by Eckern et at, the in- 
fluence of the quasiparticle dissipation is quantitatively 
weaker than that of the Ohmic dissipation in the shunt 
resistor. 10 This is due to the existence of an energy gap 
A for the quasiparticle excitation in superconductors. 
Therefore, in an ideal s-wave Josephson junction with- 
out the shunt resistance, the suppression of the MQT 
rate due to the quasiparticle dissipation is very weak at 
low temperature regime. 

In this paper, we will consider the MQT in high- 
Tc cuprate superconductor Josephson junctions. From 
many experimental studies, it is confirmed that the sym- 
metry of the pair potential (the superconducting gap) is 
d x 2_ y 2 11 ^ 12 (see Fig. 1). In such anisotropic supercon- 
ductors, the gap vanishes in certain directions (the nodal 



directions), hence quasiparticles can be excited even at 
sufficiently low temperature regime. Therefore, we will 
investigate the effect of the nodal quasiparticle dissipa- 
tion on the MQT. In the following, we will show an ana- 
lytical calculation of the MQT rate for the <i-wave c-axis 
twist Josephson junction (see Fig.l) from a microscopic 
model. Note that the effect of the quasiparticle deco- 
herence was recently discussed by Fominov et a/42 and 
Amin et al. 1A in the context of the <i-wave qubit. 

To derive an expression for the effective action, we will 
use a microscopic model of the d-wave/insulator/d-wave 
Josephson junction, described by the grand canonical 
Hamiltonian, H = Hi+H2+'Ht+'Hq, where Hi and H2 
are Hamiltonians describing the <i-wave superconductors: 
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FIG. 1: Schematic drawing of the c-axis twist Josephson junc- 
tion. 7 is the twist angle about the c-axis. Inset shows the 
pair potential for the d x 2_ y 2-w&ve superconductors. 
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where /i is the chemical potential and ip (ip^) is 
the fermion field operator. In order to obtain 
the anisotropic pair potential, the anisotropic attrac- 
tive interaction g (r — r') have to be taken into ac- 



count unlike conventional s-wave casesi* 



5 H T = 
describes 



J2 a J drdr' \t (r, r') V>L (r) yj 2a (r') + h.c. 
the tunneling of Cooper pairs and quasiparticles between 
the two sides of the junctions, and TLq = (Qi — Q2) 2 /8C 
is the charging Hamiltonian where C is the capacitance 
of the junction and Qxn) is the operator for the charge 
on the superconductor 1 (2), which can be written as 

Qi = e E<r / dr ^L ( r ) ^i<7 ( r ) ■ 

By using the functional integral methodfiSiii the 

ground partition function for the system can be written 

as follows 
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where (3 = 1/fcsT, ip(ip) is the Grassmann field which 
corresponds to the fermionic field operator ip(^), and 
the Lagrangian C is given by 

<J 1=1,2'' 

We will use the Hubbard-Stratonovich transformation 



= / VA(r,r';T)VA(r,r';T)exp 
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in order to remove the term ip in the Hamiltonian 
H(t). This introduces a complex pair potential field 
A(r,r';r) = \A(r, r'; t)| exp (icf> (r, r'; r)). The result- 
ing action is only quadratic in the Grassmann field, so 
that the functional integral over this number can readily 
be performed explicitly. The functional integral over the 
modulus of the pair potential field is taken by the saddle- 
point method. Then the partition function is reduced 
to a single functional integral over the phase difference 
<t> = 4>\ — 4>2- To second order in the tunneling matrix 
element, one finds Z = J T>4>(t) exp [— <S G ff /ft] , where 
the effective action is given by 
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This expression coincides with the previous resultsiiSiiS 
The second term in eq. © describes dissipation due to 
the quasiparticle tunneling. The third term describes 
the tunneling of Cooper pairs (the Josephson tunneling) . 
The kernels a(r) and /3(t) are given in terms of the diago- 
nal and off-diagonal components of the Matsubara Green 
function in Nambu space, usually denoted by Q and J- 

Q ( T ) = -f £I*(M')| 2 ft (*:',-•>-), (6) 

k,k' 

^ = -|E |^fc')| 2 ^i(fc^)^2 t (fc',-T). (7) 
k,k' 

The Green functions are given by 



Q (k,u) n ) = 



ft (ihw n + Cfe) 
~(^„) 2 +£| + A(fc) 2 

ftA(fe) 
(ftu;„) 2 +£ 2 +A(fc) 2 ' 
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where £u — h 2 k 2 /2m — fi and ftw„ = (2n + 1)tt/P is 
the fermionic Matsubara frequency (n is an integer) . In- 
formation about the anisotropy of the pair potential is 
included in A(fe). In the case of the cuprate high-T c su- 
perconductors (the d x 2_ y 2 symmetry), A(fe) = Aocos20 
(see Fig .1). 

We now turn to the calculation of the effective action 
and the MQT rate (the inverse of the lifetime of the 
metastable state) for the c-axis twist Josephson junction. 
In Fig.l, we show schematic of this junction. In this fig- 
ure, 7 is the twist angle about the c-axis (0 < 7 < tt/4). 
Such a junction was recently fabricated by using the 
single crystal whisker of Bi2212 (BiaSraCaC^Os+s). 20 
Takano et al. measured the twist angle dependence of the 
c-axis Josephson critical current and showed a clear evi- 
dence of the d x 2_ y 2 symmetry of the pair potentiali22*2i 

In the following, we assume that the tunneling between 
the two superconductors is described in terms of the co- 
herent tunneling (|i(fc,fe')| = \t\ Sj^^ ^ , where fe|| is 

the momentum parallel to the a6-plane.). For simplic- 
ity, we also assume that each superconductor consists 
of single Cu02 layer, Ai(fc) = Aocos26', and A2(fe) = 
Ao cos 2 (9 + 7). Moreover, we consider the low tempera- 
ture limit (fcsT <C Ao). In the case where the Josephson 
junction is biased by an externally applied current l e xt, 
we have to add an additional potential contribution linear 
in 4>m^ At this level of approximation, the effective action 
of the current biased c-axis twist Josephson junction is 
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where M = C(fi/2ef 
washboard potential 



is the mass and U{<j>) is the tilted where 



U{<t>) = -Ej(j) cos0 + 



(12) 



In this equation, Ej = (H/2e) Iq is the Josephson cou- 
pling energy and Ic = —(2e/H) J Q h/S dr/3(r) is the Joseph- 
son critical current. From eqs.(7) and (9), we can obtain 



lc(0) = f |*| 2 iV Q 2 Ao, 
Zc(7r/8) « 0.66J c (0), 



(13) 
(14) 



where iVo is the density of states at the Fermi energy. 
The result of /c(0) agrees with that of Tanaka and Kashi- 

29 

wayaifi* 

In the following, we will consider the effect of the nodal 
quasiparticles on the MQT. For this purpose, we first 
calculate the dissipation kernel a(r) for two types of the 
c-axis junction, i.e., (1)7 = and (2) 7 7^ (here we will 
show the result for 7 = tt/8 case only.). Note that the 
behavior of a(r) have been already predicted in Ref. [18] . 
Below we will derive the analytic expression for a(r) and 
calculate the renormalization mass. 

In the case of the c-axis junction with 7 = 0, the nodes 
of the pair potential in the two superconductors are in the 
same direction. Therefore, the node-to-node quasiparti- 
cle tunneling is possible even at very low temperatures. 
In this case, the asymptotic form of the dissipation kernel 
at the zero temperature is given by 



a{r) 



3ft 2 |<| 2 iV 2 1 
16A |r| ; 



(15) 



for Ao|T|/ft ^> 1. This gives the super-Ohmic dissipation 
which agrees with Ref. [18] and [19]. Note that the dis- 
sipation kernel for the normal Ohmic shunt resistance is 
a (r)~l/r 2 i2^ 

On the other hand, in the case of the finite twist angle 
(7 = 7r/8), the asymptotic behavior of the dissipation 
kernel is given by an exponential function due to the 
suppression of the node-to-node quasiparticle tunneling, 



i.e., 



a(r) ~ exp 



1 AqM 
V2 n 



(16) 



for Aq\t\/H ^> 1. This result coincides with the previous 
prediction^ Eq. (16) is very similar to that of the con- 
ventional s-wave Josephson junctions with the constant 
pair potential A: a(r) ~ exp (— 2A|r|/7i) . — If the phase 
4>{t) varies slowly with time on the time scale given by 
H/A.Q, we may expand 0(t) — </>(t') in eq.(ll) about t—t'. 
This gives 



SM(ir/8) 



(17) 



6M(n/8) fa 0.25 



2 N§ 
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(18) 



Hence under above condition, the dissipation action 
acts as a kinetic term so that in the end the effect of 
the quasiparticle dissipation results in an increase of the 
mass, i.e., M — > M + SM (the mass renormalization). 

The MQT rate at the zero temperature is given by the 
formulai* 2 ^ 



r = lim — ImlnZ. 



(19) 



In order to determine F we employ the instanton the- 
oryi 2 ^ When the barrier is low enough for the MQT to 
occur but still so high that the instanton approximation 
is valid, Y is given by 



^4exp — 



Sr- 



(20) 



where Sb — 5eff[<As] is the bounce exponent, that is 
the value of the the action S e g evaluated along the 
bounce trajectory 4>b(t). Using above and Leggett 
et.al.'s methodjSi we obtain the main results of this pa- 
per, namely, analytic expressions for the MQT rate (Note 
that we have set 7 e xt ~ Ic so that we can approximate 
the washboard potential U(<ft) as a quadratic-puls-cubic 
one.): 
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(21) 



(22) 



where, 



B( 7 ) 



1 - 



led) 




(23) 



and To (7) is the decay rate without the quasiparticle dis- 
sipation. In eq.(23), 5M(n/8) is given by eq. (18) and 



^ 2 A 



vT 



X 2 Jq 



dss 1 K 1 {sx) 2 {2A) 



where uj p is the plasma frequency and K\ is the modified 
Bessel function. As an example, for An = 42.0 meV, 25 
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7 C ( 7 = 0) = 1.45 x 10- 4 A, C = 10 £F, and I ext /Ic{l) = 
0.9, we obtain 

£(7)_ f 90% for 7 = 

r ( 7 ) ~ \ 96 % for 7 = tt/8 ' lZ0J 

Therefore, the node-to-node quasiparticle tunneling in 
the case of the 7 = junction gives rise to large reduction 
of the MQT rate in compared with the 7 = tt/8 case. 

To summarize, we have presented the analytical cal- 
culation of the MQT rate for the c-axis twist Josephson 
junction by making use of the functional integral method 
and the instanton theory. Within the coherent tunneling 
approximation, we find the super-Ohmic dissipation in 
the case of the zero twist angle junction. This dissipa- 
tion is caused by the node-to-node quasiparticle tunnel- 
ing between the two superconductors. In the case of the 
finite twist angle, on the other hand, the suppression of 
the MQT rate is weak in compared with the 7 = case 
due to the inhibition of the node-to-node quasiparticle 



tunneling. 

In this paper, we have considered the c-axis Joseph- 
son junctions. In <i-wave Josephson junctions along the 
afc-plane (e.g., YBCO/PBCO/YBCO ramp-edge junc- 
tions^ and YBCO grain boundary junctions^), the zero 
energy bound states (ZES) 28 29-30 give a crucial contri- 
bution to the Josephson and the quasiparticle current. 
Therefore, it is very interesting to investigate the effect 
of ZES on the MQT from a microscopic Hamiltonian. 

Finally we would like to point that, in a phase type 
qubitj 3 ^* 3 ^ the MQT is used in final measurement process. 
We expect that our result will help in understanding the 
decoherence in the measurement processes for the d-wave 
phase qubits. 
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